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PROJECTED HARTREE-FOCK SPECTRA WITH 


WOOD -SAXON BASIS FUNCTIONS 
by Richard C. Braley and William F. Ford 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

We have previously reported the results of an investigation of intrinsic 
Hartree-Fock spectra and projected nuclear properties, such as E2 rates 
and inelastic scattering cross sections, when basis functions with Wood- 
Saxon radial dependence are used. It was found that predicted values of 
properties sensitive to the long range behavior of the wave functions were 
improved by the use of a Wood-Saxon basis. Since the basis functions are 
normalized to unity, the short-range behavior of the wave function, and 
hence the energies, should be affected. This point is investigated by com- 
parison of projected Hartree-Fock spectra using basis functions having 
harmonic oscillator and Wood-Saxon radial dependence. 

INTRODUCTION 

The study of nuclear properties based on the Hartree-Fock method 
has been limited almost entirely to calculations made with harmonic oscil- 
lator (HO) model spaces. This is true in spite of the well-known incorrect 
asymptotic behavior cf the basis states. Recently, the present authors in- 
vestigated the use of Wood-Saxon (WS) basis functions in Hartree-Fock 
studies. ^ It was found that the use of Wood-Saxon functions yielded im- 
proved results for some observables relative to those obtained with HO 
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functions. ’ Those calculations were made with an inert O core, and 

only the 2s -Id shell model states were used as the basis space. It has 
been found, however, that core polarization effects make important con- 
tributions to many nuclear properties when HO functions are used and thus 
such effects should be investigated for the WS case also. In the earlier 
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work the Wood -Saxon results were compared with similar ones obtained 
using HO functions. The choice of oscillator parameter posed something 
of a problem but at length two physically meaningful parameters were 
chosen. One of the parameters (b = 1. 69 fm) yields peaks and nodes 
which match those of the Wood -Saxon functions; this is probably the most 
logical choice for purposes of compcrison. The other (b = 1.924 fm) was 
chosen to reproduce approximately binding energies, gaps, and intrinsic 
single particle energies obtained in the WS results. 

Since the WS functions do not fall off as rapidly as the HO functions, 
one would expect to find an improvement in the nuclear properties which 

are sensitive to the long range behavior of the nuclear wave function. 
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This was found to be the case for the study using an inert core. However, 
the long range behavior should also manifest itself in properties such as 
energies which are sensitive to the short range behavior of the wave func- 
tion. This is due to the fact that the single particle basis functions are 
normalized to unity. 

In this article we compare the results of HF calculations made with 
the three model spaces defined above. The assumption of an ^ O core 
is dropped, making it possible to study the effects of core polarization. 
Several projected properties such as RMS radius, E2 rates, and projected 
energy spectra are examined and compared for the three spaces in an 
effort to determine any effect which long range properties may induce in 
the energy spectra. 

Calculation 

The HF model space used in the calculations consisted of the shell 

model states: ^Pi/2» ^^3/2* ^^ 5 / 2 * ^^3/2' 

evaluation of the 2 -body matrix elements for the harmonic oscillator 
functions is straightforward; however, when WS functions are used, 
their evaluation requires special techniques which have been discussed 
in an earlier publication. The nucleon -nucleon interaction used in these 
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calculations is the Volkov force 
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V(|?j-r2|)= {a+ -T2 + 6ctj ■ Vj • TjjfClFj - TjI ) 


where 


£(|Fi-J"2 


= -78 e 


-dv?2l^ 


/2. 25) 


+ 82. 5 e 


-(|rj-r2l^/0.64) 


This force (with a = 0. 1875, ^ = 0. 0775, y = 0. 2025, 6 = 0. 1775) has 
been found to yield good saturation properties for the p- shell nuclei. 

Most of the nuclear properties which we shall be concerned with in- 
volve the evaluation of reduced matrix elements (between projected HF 
states, sums of single-particle operators which may 

be expanded in terms of spin-angle tensors with coefficients i^L 3 j(r). 
Reduced matrix elements of such an operator may be expressed as^ 


•III wIK)= 2 f 

n=l / LS JO 


OO 


.if 


(r)w. 




r(r)r^ dr 


where the structure information is contained in the transition densities 
p^Sj(r). Of course the calculation of the projected energy spectrum also 
involves two-body operators, and is therefore more complicated. One 
must calculate 


where the Hamiltonian, H, is the sum of one-body (kinetic energy) and 
two- body (potential) operators. The calculation of this quantity is dis- 
cussed in detail elsewhere. ^ 


4 


RESULTS AND DISCISSION 

The radial functions for sonr*e of the single- particle states are com- 
pared in Figures 1-3. As noted earlier the WS functions are generally 
depressed (relative to the HO) it small values of r, but fall off less 
rapidly at larger values. This is a consequei'ce of normalization to 
unity. 

2T 24 28 

The intrinsic and projected quantites for Ne, Mg, Si, and 
32 

S are presented in Tables I and n. First of all vve note that in all 
cases the HO results yield a lower HF energy (more binding) than the 
WS. This is due to che fact that the short range part of the WS function 
is depressed relative to that of the HO type. In Table I, we note a sig- 
nificant difference in the intrinsic properties (Ejjp and A), as was 
found for the inert core study. A similar difference is noted for the 
projected properties, but comparison with Table HI reveals that the WS 
rpsuits are in better agreement with experiment. In Table H (b = 1. 924 fm. ) 
the intrinsic and projected properties are seen to be in much closer agree- 
ment than in the previous case. However, the radii are generally too 
large for this oscillator parameter. 

The projected spectra for the nuclei of interest are presented in 
Figures 4-7. In all cases, the spectra for the two oscillator parameters 
differ very little; furthermore, the Wood-Sa . spectra are only slightly 
different from the HO results. Comparison with experiment reveals that 
the predicted spectra in all cases are too compressed. This is not sur- 
prising, however, due to the small size of the model space. 

20 

The Ne Wood-Saxon result show a somewhat peculiar distribution 
of the different amounts of the various J-components in the intrinsic state. 
This gives rise to an unusual spectrum, which does not appear to be ro- 
tational. The wave function for this state v/as examined very closely for 
errors and appears to be correct. It may be that the particular solution 
obtained is not the best HF self-consistent solution >vhich can be found, 
but represents only a local minimum of the HF energy surface. 
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CONCLUDING REMARKS 

The purpose of this work has been to determine whether the use of 
WS basis functions in HF calculations would affect the projected energy 
spectra. Previous studies have established that the long range behavior 
of these functions modified the single particle spectra slightly and lowers 
the HF energy relative to HO results. Our present findings show that 
absolute values of the energy spectra are higher for the WS functions but 
that the relative) separations of the projected energy levels are not affected 
significantly. 

However, this study has provided an additional valuable result, namely 
that one need not necessarily examine such hard- to- calculate properties 
as projected energies in order to test the relative validity of the various 
aspects of a nuclear model. Due to the small number of basis states 
in the model space it was not expected that the projected spectra would 
agree very well with experiment. However, it was rather surprising 
to find that the large differences in intrinsic states and single particle 
properties did not manifest themselves significantly in the energy spectra. 

It appears that an efficient procedure to follow when the standard HF ap- 
proach is used is first to examine the effects which model changes have on 
single particle properties (radii, E2 rates, inelastic nucleon scattering, 
etc. ), and then when one is satisfied that the moder wave functions have 
improved significantly, to examine properties such as projected spectra 
which are much more involved in terms of the calculational tools required. 
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TABLE I. - OSCILLATOR PARAMETER, b = 1.69 fm 


Nucleus 

Basis 


20xro 

HO 

2.65 

Ne 

WS 

2.84 


HO 

2.75 

WS 

2.96 

28 si 

HO 

2.81 

WS 

3.09 

32g 

HO 

2.86 

WS 

3.17 


B(E2; O'^ - 2 *), 
. fm^ 

MeV 


77.3 

-118.43 

10.6 

140.8 

-84.98 

8.0 

111.2 

-139.31 

0.38 

174.4 

-87.54 

0.11 

157.6 

-192.15 

11.51 

274.9 

-116.03 

8.6 

112.4 

-232.83 

0.38 

182.7 

-134.15 

0.22 


TABLE n. - OSCILLATOR PARAMETER, b = 1.924 fm 


Nucleus 

Basis 


B(E2; 0*^ - 2'^), 
e2 . 

®HF» 

MeV 

n 


HO 

3.01 

135.7 

-98.37 

8.3 

WS 

2.84 

140.8 

-84.98 

8.0 


HO 

3.12 

191.2 

-116.77 

0.24 

WS 

2.96 

174.4 

-87.54 

0.11 

00 

CO 

HO 

3.20 

269.8 

-160.15 

9.1 

WS 

3.09 

274.9 

-116.03 

8.6 

32g 

HO 

3.25 

192.8 

-193.75 

0.24 

WS 

3.17 

182. 7 

-134.15 

m 




































































TABLE in. - EXPERIMENTAL 


E2 RATES AND RADn 
























